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ABSTRACT. Medium-resolution (8v=45 km s' 1 ) optical spectra of the bright, high-redshift (z=4.1) quasar 
0000-263 taken at the ESO 3.5-m NTT telescope were analyzed to determine the distribution of column densities, 
velocities and line widths of the Lyman-a forest absorption components. The values of Nj-j, b and z were determined by 

fitting Voigt profiles to the lines, and convolving with a Gaussian instrumental response function. Over 350 components 
with log Nh > 13.2 were identified. An analysis of the dependence of the number of components with z reveals that the 

number evolution of components obeys the power law dN/dz « (l+z)Y , where y =0.5 ± 0.4 for the sample of 1 82 lines 
with log Npj > 14.0. The distribution of component strengths is found to obey f(Njq) Npj'P, where P = -1.55 for 
components with log(Np[) > 14.7, and P = -0.68 for the components with log(Npj) > 13.5. A distinct break in the f(Nj-{) 
histogram is also observed, at log(Np[) - 14.7. The results are briefly considered in the context of theoretical models of 
quasar Lyman alpha clouds and their evolution. 

OBSERVATIONS. Optical spectroscopic observations of the quasar 0000-263 were obtained during two nights 
at the ESO 3.5-m NTT telescope at Cerro La Silla, Chile. Coadded CCD images of the echelle spectra were flat fielded, 
cosmic rays were removed, and spectra were extracted and wavelength calibrated. The average FWHM resolution of the 
spectrograph was 45 km s*l, and the spectrum was sampled at 0.3 A Intervals. The S/N for the coadded spectrum is 
estimated to be 10, based on analysis of the line-free region redward of the QSO emission line at 6200 A. Line-free points 
of the spectrum blueward of 6200 A were selected and a 4th degree spline was fitted to remove the continuum. The 
spectrum, and the continuum fit used are shown in Figure 1 . The Lyman-a forest was then analyzed by fitting Voigt 
profiles convolved with the instrumental response of the spectrograph to the absorption components. The fit chosen was 
based on the minimization of the residual between observed data and modelled profile, after iterating the parameters Nj^.b, 

and v for each component from an initial component model which was determined interactively. 356 components were 
found, of which 182 had log(N) > 14.0. The sample should be complete for components with log(N) > 14.0, which is a 
conservative completeness limit corresponding to an equivalent width threshold of - 700 mA. The components were 

analyzed to determine if there were statistically significant correlations in the column densities N and Doppler parameter b, 
and to determine evolutionary trends in dN/dz and in the distribution function of component strengths, f(Nj^). 

RESULTS. The correlation between N and b values yields a correlation coefficient of R=0.47. The correlation is 
statistically significant but may be affected by systematic effects intrinsic to line profile modelling of spectra which have 
many blended components, such as those described by Rauch, et al (1992). After excluding components with < 0.65 

A, and b > 90 km s'l, the correlation coefficient was found to decrease to R=0.25. This suggests that the excluded 
components were responsible for much of the appearance of correlation, although a correlation between b and N for weak 
components has been predicted from theoretical models of differential ionization of Lyman alpha clouds (Donahue and 


Figure 1: The coadded spectrum of QSO 0000- 
263, with the continuum fit superimposed. Lyman 
a components were fit from the wavelengths 
5400A to 6200A, using modelled Voigt profiles 
which were convolved with a Gaussian 
instrumental resolution function, which had a 
FWHM of 5v=45 km s'l. 


5000 5500 6000 6500 7000 7500 

Wavelength 

377 


Shull, 1992). 
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The distribution function of the component column densities f(Nfj) was examined, and the values of f(Npj) for 
log(N H ) > 13.75 are presented in Figure 2. The function f(N^) is usually described in terms of a power law f(Npj) « 
NpfP, and for the current work we find a value of P = 1.55 for the strongest components with log(Npj) > 14.7, and a 
lesser slope of p = 0.7 for the entire sample. The value of P reported here for the stronger components is consistent with 
other Lyman a systems, where P is typically 1.5-1.7, (Charlton, et al 1992). While the break in the power law has been 
predicted by theoretical models to result from changing ionization state of the Lyman a clouds, it may also be due to the 
weaker lines being incompletely detected due to S/N limitations. 

To study the evolution of the number of Hydrogen components N com p S with redshift z, we considered redshift 
intervals of Sz = 0.05 and examined the number of components in each interval. A plot of N vs. (1+z) is shown in Figure 
3, for the components with log (Npj) > 14.0. A power law fit was applied to the data, and it was found that dN/dz « 
( 1 +z)Y, where y =0.5 ± 0.4. The fit is shown on Figure 3 with the dashed line. The value of y derived here is much less 
than the average value of y = 2.36 ± 0.4 from Bajtlik, Duncan and Ostriker (1988), based on several quasar spectra at the 
lower redshift redshift range 3.7 > z > 1.7. The value of y =0.5 ± 0.4 is however consistent with Lyman limit systems, 
for which y =1.1±0.4 for systems with redshifts ranging from 4 > z > 0.4 (Charlton, et al 1992). It is also worthwhile to 
note that the value of y =0.5 is consistent with the cosmological model of a constant number of components per unit 
comoving volume (Sargent, et al 1980; Charlton, et al 1992) where dN/dz « (\+z)I(\+Hqz)^ 1 ^-. For Oq= 1, and assuming 
pressure and ionization rates constant over the redshift interval considered, the model reduces to dN/dz °c (l+z)0-5 ( 
consistent with our observed result. Further observations at higher S/N and at higher spectral resolution are needed to 
determine if the apparent low value of y for QSO 0000-263 is indeed real, and whether it does indicate a change in the 
number evolution of Lyman alpha clouds at high redshift. 
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Figure 2: (left) Plot of the histogram function log(f(NH)) with power law fits f(Njj) Npj'P, with b=1.55 and 0.66. 
Figure 3: (right) Plot of the number of components N^^ps against (1+z), with the best fit dN^j^pg/dz « (l+z)0-5 
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